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Out-of-equilibrium Anderson model at high and low bias voltages 

Akira OGURI 

Department of Material Science, Osaka City University, Sumiyoshi-ku, Osaka 558-8585, Japan 

(Received May 22, 2002) 

We study the high- and low-voltage properties of the out-of-equilibrium Anderson model 
for quantum dots, using a functional method in the Keldysh formalism. The Green's function 
at the impurity site can be regarded as a functional of a nonequilibrium distribution function 
f c g(u>). The dependence of the Green's function on the bias voltage V and temperature T arises 
through /efi(w). From this behavior as a functional, it is shown that the nonequilibrium Green's 
function at eV — > oo is identical to the equilibrium one at T — > oo. This correspondence holds 
when the couplings of the dot and two leads, at the left and right, are equal. In the opposite 
limit, for small eV, the low-energy behavior of the Green's function can be described by the 
local Fermi-liquid theory up to terms of order (eV) 2 . These results imply that the correlation 
effects due to the Coulomb interaction U can be treated adiabatically in the two limits, at high 
and low bias voltages. 
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1. Introduction 

The Kondo effect in quantum dots is a very active field 
of current research, and recent experiments 1-4 ) have been 

, shown to be in qualitative agreements with early pre- 
dictions. 5 ~ 7 ) Theoretically, the equilibrium and linear- 
response properties of the Kondo system realized in the 

' quantum dots have been well understood through those 
of magnetic alloys 8 ) and precise calculations with the nu- 
merical renormalization group method. 9 ) However, the 

, nonequilibrium properties under a finite bias voltage V 

• have not yet been fully understood. It is a novel problem 
of the strongly correlated electron systems, and has been 
studied extensively. 10 ~ 14 ) 

At small voltages eV <C Tk near the linear-response 

■ regime, nonequilibrium properties at low energies, i.e., 
at T -C Tk and uj -C Tk, can be described by the local 
Fermi-liquid theory, 15 ' 16 ) where Tr is the Kondo tem- 

, perature. Specifically, the nonlinear-response of the cur- 

■ rent through the quantum dots has been calculated up 
to terms of order V 3 based on the Kondo model 13 ) and 
Anderson model. 14 ' The coefficients can be expressed 
in terms the correlation functions defined with respect 
to the equilibrium ground state, 14 ' and in the electron- 

, hole symmetric case the coefficients are determined by 
Tk and the Wilson ratio R. In this paper, we suggest 
a procedure to estimate Tk and R experimentally from 
the differential conductance near the unitarity limit. 

The main purpose of this paper is to describe a rela- 
tion between the effects of the bias voltage and temper- 
ature based on a property of the generating functional 
Zj for the Keldysh Green's function of the Anderson 
impurity. We show that Zj can be regarded as a func- 
tional of a nonequilibrium distribution function f e g(ui), 
and the dependence of Zj on eV and T arises through 
/ e ff(aj) . From this feature of the generating functional, it 
is deduced that the nonequilibrium Green's function at 
eV — ► oo is identical to the equilibrium one at T — ► oo. 
This correspondence holds exactly when the couplings 
between the dot and two leads are equal. Furthermore, 
the low-voltage Fermi-liquid behavior mentioned above 



can also be deduced from this property of Zj. As the 
functional of Luttinger and Ward 17 ) has played an im- 
portant role in the usual Fermi-liquid theory, some of the 
nonequilibrium properties can be deduced from those of 
the generating functional. 

Our results obtained in the two opposite limits of eV 
imply that the Coulomb interaction U can be treated adi- 
abatically at both low (eV < T K ) and high (U < eV) 
bias voltages. Therefore, as one of the possibilities, we 
could expect that the perturbation theory in U works 
for all values of eV, although another possibility that a 
different phase of non-perturbative nature exists at in- 
termediate values of eV could not be ruled out. 

2. Keldysh Formalism for Anderson Model 

2.1 Green's junctions 

We start with the single Anderson impurity connected 
to two reservoirs at the left (L) and right (i?): 

H= Hi + H 2 , (1) 

\=L,R ko a 

H 2= XS (4c Act + 4x) + Urmnu , (3) 

X—L,R a 

where d a annihilates an electron with spin a at the dot, 
and nda = d\d a . In the leads A (= L, R), the excitations 
are described by = tk + eV\ and corresponding eigen- 
function 4>k\( r )- To specify the static potentials V\, we 
introduce additional parameters ct\ as Vl — c^lV and 
Vr = —olrV with <xl + otR — 1- Here, the Fermi level at 
equilibrium is taken to be the origin of the energy. The 
onsite potential Ed is assumed to be a constant which 
does not depend of the bias voltage. The mixing matrix 
elements v\ describe the couplings between the dot and 
leads, and c A(T = c fcAcr </> fcA (r A ) annihilates an electron 
at the interface r x . We assume the local density of states 

l 1 2 

PaM = Efe |0feA C^a) I S (v - e k\) to be a constant and 
its band width D to be infinity. We will use units h = 1 . 
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Fig. 1. The Keldysh contour of the time evolution. 



In order to describe a nonequilibrium steady state re- 
alized under the applied bias voltage, we employ the 
Keldysh formalism. 18 ) We define four types of local 
Green's functions at the impurity site; 



G—(t) = -i<T (*„(*) 4(0)) , 
G-+(t) = i(4(0R(t)>, 
G + -(t) = -i(d a (t)4(0)), 
G++(i) = -i<fd„(t)4(0)>, 



(4) 
(5) 
(6) 
(7) 



where T and T denote the time-ordering and anti-time- 
ordering operations, respectively, and d a (t) is a Heisen- 
berg operator. Note that these functions are linearly 
dependent G"+ + G+" = G — + G++, and the re- 
tarded and advanced Green's functions are given by 
G r = G— - G~+ and G a = G~ - G+~. The average 
(• • • ) is taken over the density matrix p(t) at t = 0. The 
initial condition is given at t — — 00, where the two leads 
are separated from the dot and each part is in its own 
thermal equilibrium described by H\ and the chemical 
potentials /i\ = eV\. Then H2, which includes the mix- 
ing v x and interaction U, is switched on adiabatically. 
The time evolution can be treated in the interaction rep- 
resentation (see Appendix A), where the Wick's theorem 
is applicable for the time-ordered correlation functions 
along the Keldysh contour shown in Fig. 1 . 

In the noninteracting case U = 0, the Green's func- 
tions which include all contributions of the mixing v x 
are given, as functions of the frequency w (see Appendix 
B),by 

Go - M - [1 - ./cffM] G5M + jcffM G » , (8) 

G o + ( u> ) = ~ /cff M [ G5 M - Gg M ] , (9) 
G+-(oj) - [1 - /cffM] [GSH - Gg(w)], (10) 

G++{w) = - [1 - ./effH] GS( W ) - /effM W • (11) 

Here G|J(w) = [w-^ + iA]" 1 and G a (uj) = {G r {u)}*. 
The level width is A = T L + T R with T\ = 7rp A vf . The 
function f c e(uj) expresses the nonequilibrium distribu- 
tion of the electrons at the impurity site, 11 ) 



/eff(w) = 



fL(u)T L + f R (w)r R 



(12) 



where f x (w) = f{u - fi x ) with f(u) = [e u ' T + l]" 1 . At 
T = 0, f e s(u>) has two steps at o> = (1l and as shown 
in Fig. 2. Note that the Keldysh formalism is applicable 
also in equilibrium, i.e., at eV — 0, where / e ff(w) becomes 
equal to the usual Fermi function f(u>). 

The Green's function for U 7^ satisfies the matrix 



/effM 



Fig. 2. Nonequilibrium distribution f e s(u>) at T = 0. 

Dyson equation; {G(w)} _1 = {Go(w)} _1 — E(u) , 

"Gn 



Go = 



£ = 



s— 



(13) 



Here S(w) is the self-energy due to the interaction U. 
Four types of the self-energies are also linearly dependent 
£"+ + S+- = -S — - £++. The perturbation theory 
for G(u>) is described in the real-frequency (or real-time) 
representation, and thus the dependence of G(w) on eV 
and T arises only through f e s(ui) which enters the non- 
interacting one Go(w). This feature seen in the Keldysh 
formalism is quite different from that of the Matsubara 
formalism in which the temperature dependence arises 
through the summations over the imaginary frequencies. 

2.2 Generating Functional 

In order to see this feature of the Keldysh formalism 
more explicitly, we employ the generating functional Zj 
that yields the perturbation of series for G(u). It is 
given, in the path integral form, 19 ) by 

Zj = JvrfVv e is ^>n) ; ( 14 ) 

Sfat.rj) = S (v\v) + Sex^r?) + Suiv^rj). (15) 

The action S consists of three parts corresponding to the 
free, external-source, and interaction contributions; 

/oo 
dtdt' V Ut)K (t,t') Vr7 (t'), (16) 



SuW,v) 



t ^ - 



/oo 
-oc 



dt 



v Ut)JAt) + Ji(t)vAt') 

(17) 

rj\_(t) v ,_(t) v\.(t) Vl _(t) 



v\+{t) Vu (t) V \ + {t) Vl+ {t)]. (18) 



Here r)\ (t) = ^ 77* _ (t) , n\+ (t) ^ is a two component field 

of the Grassmann number. The label =F specifies the 
branches of the Keldysh contour (see Fig. 1) where each 
of the components corresponds to. The Kernel K in eq. 
(16) is the inverse matrix of the noninteracting Green's 
function, 



K (lj) = {Go(^)}- 1 

dui 
2^ 



K (t,t') = / 

J — < 



(19) 

Koi^e-^-t'K (20) 
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In eq. (17), jt(t) = is an external 

source of the anticommutating c-number, which is intro- 
duced for the later convenience. In eq. (18), the sign of 
the first and second terms are determined by that of the 
exponent of the time-evolution operators £/(+oo,— oo) 
and U(— oo,+oo), respectively (see Appendix A). For 
U = 0, the path integral of the fermionic fields in eq. 
(14) can be evaluated analytically 



VrfVrj e 



i[So(.) t ,»))+S cx (r ) t >77 )] 



z j=o exp 



/OO 
dtdt'4(t) GoMJ.tf) 
-OO 



(21) 



where Zj =0 = limj^o-^j- Then following the stan- 
dard prescription, 19 ^ the generating functional and full 
Green's function can be rewritten in the form 



Z, 



7° 



Zj SJL(t) 5J ^(t') 



z., 



(22) 



, (23) 



.7=0 



where v, v' = +, — . In eq. (22), the fermionic fields 
in Su have been replaced with the functional differen- 
tiations, ry av => —iS/dJcrv and r\ av => iS/SJ^. The 
perturbation series of G can be generated from eq. (23) 
by expanding e lSu in eq. (22) as a power series of Su 
and then carrying out the functional differentiations us- 
ing the explicit expression of Zj given by eq. (21). Each 
term of the perturbation series is written in terms of Go 
and the integrations over the internal variables of real 
time (or real frequency). Therefore, Zj can be regarded 
as a functional of / c g(w), and the dependence of G(w) on 
eV and T arises only through f c s(uj). From this of prop- 
erty, the behavior of G(oj) at high and low bias voltages 
can be deduced exactly. 

3. High- Voltage Behavior 

3. 1 General properties 

We consider the behavior of G(ui) in the high-voltage 
limit eV — > oo. To be specific, we assume that the bias 
voltage is applied to be oil > and an > 0. In the limit 
of eV — > oo, the chemical potentials tend to (il — » oo 
and fiR — > — oo, and thus the distribution function be- 
comes a constant fesiu)]^^^ = T L /(T L +T R ), see 
Fig. 2. When the couplings with the two leads are equal 
= Vr, the value of this constant becomes 1/2 and 
/ e ff(<*')| e v- > cx> coincides with the high-temperature limit 
of the usual Fermi function /(w)| r _, = 1/2. Therefore, 
the noninteracting Green's function Go for eV — ► oo be- 
comes identical to that for T — ► oo at equilibrium. The 
same correspondence holds also for the full Green's func- 
tion G since it is determined by eqs. (21)-(23) for given 
Go- Thus in the case of Tl = Tr the two limits, i) 
eV — » oo and ii) T — > oo keeping eV = 0, are equiva- 
lent as far as the Green's function at the impurity site 
is concerned. Note that the temperature is not neces- 
sary to be kept at T = in the eV — ► oo limit. In 
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Fig. 3. Spectral function obtained with order U 2 self-energy. 



these two equivalent limits, the Coulomb interaction U 
may be treated adiabatically. This is because in thermal 
equilibrium the perturbation series in U is absolutely 
convergent at T = for any finite J/, 20 -* and there is 
no phase transition at finite T. When the voltage is fi- 
nite but still much larger than the other energy scales 
eV > max({7, \E d \, A,T), the behavior of G(u) at low 
frequencies <C eV may be approximated reasonably 
by that in the eV — > oo limit. 

In the case of Tl =/= the two limits, i) and ii), are 
not equivalent because the charge distribution around 
the impurity at eV — > oo does not correspond to that 
at T — > oo. In this relation, we examine the T — > oo 
limit at finite eV. In this limit the distribution function 
is given by / c fr('-^)lT^oo = 1/2 without the assumption 
of Tl = Tr, and the Green's function becomes identical 
to that in the T — > oo limit at eV = 0. Physically, this 
is rather obvious because at T ^> eV the effects of the 
bias voltage are hidden by the thermal fluctuations. 

3.2 Order U 2 self-energy 

In order to show a rough sketch of the Green's func- 
tion in the high- voltage limit, we use here the order U 2 
retarded self-energy. 11, 16, 21 -* It yields a reliable picture 
at least qualitatively, and can be calculated analytically 
at eV — * oo, 



i 



i3A 



(24) 



Here we have assumed the electron-hole symmetry tak- 
ing the parameters to be = —U/2, [i^ = eV/2, 
fj,R = —eV/2, and Tl — Tr (= A/2), so that the results 
are applicable equally to the T — > oo limit at equilibrium. 
Note that eq. (24) is the correlation part, which is ob- 
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tained by separating the Hartree-Fock contribution U/2 
as E r = U/2 + X£ or , and the retarded Green's function 
is given by G t (uj) = [u + iA — EJ^u;)]- 1 . Fig. 3 shows 
the spectral function, ImG"*(o;), obtained using the order 
U 2 self-energy. The Kondo resonance situated at co = 0, 
eV = (dashed line) disappears at high voltages (solid 
line), as seen in Fig. 3(a) for U/(wA) = 2. The results 
in the eV — > oo limit are plotted for several values of U 
in Fig. 3(b). For strong interactions U ^> A the spectral 
function has two peaks corresponding to the Hubbard 
bands at ui ~ ±{7/2, while for weak interactions U < A 
the contributions of the mixing dominate and it results 
in the single-peak structure. The order U 2 results could 
be refined quantitatively by including the higher-order 
terms 22 ) or by numerical methods. 
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Fig. 4. The U dependence of R and A : Bethe ansatz results. 



4. Low- Voltage Behavior 

4-1 Ward identity 

In the opposite limit, at small bias voltages eV <C Tk, 
the low-energy properties can be described by the Fermi- 
liquid theory. 15 ' 16 ) The proof has been provided in 
the previous paper using the Ward identity for the first 
and second derivatives of the self-energy with respect to 
e y. 14 ) We describe here the outline briefly to emphasize 
the properties of S(w) as a functional of the distribu- 
tion function f a s(u>). Since the voltage V enters Gq(u>) 
through /eff(w) as seen in eqs. (8)-(12), we have 
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d(eV) 

d 2 
d(eV) 



v=o 



d_ 

du> 



d_ 

dui 



d 
dE~d 

d 



8E d 



G0:eq(w), (25) 



G 0:eq (^). (26) 



Here Go :eq (w) = Gq{uj)\ v=0 is the equilibrium Green's 
function, a = (q^lTl — cbJ"r)/(^l + Fr), and k = 
(a|r L +a 2 R T R )/(T L +T R ). The derivatives of S(w) 
with respect to eV can be calculated by taking the 
derivative of Go in each term of the perturbation se- 
ries, and then replacing the derivative d/d(eV) with 
(d/duj + d/dE d ) using eqs. (25) and (26). It yields 



d(eV) 
5 2 S(w) 



d{eVf 



v=o 



v=o 



d_ 

d_ 



d 

d 
dE~d 



S eq (w) , (27) 



Seq(w) 



7J 2 S cq (w) 



(28) 



where £ cq (w) = S(w)|y =0 . In eq. (28), D 2 denotes the 
functional operation carrying out the second derivative 
(d/du> + d/dEd) 2 for all the single Go's in the perturba- 
tion series of S cq , which can formally be expressed using 
the functional differentiation of S cq with respect to Go, 



d d 
du/ + 'dE„ 



GZAu') . (29) 



Note that the perturbation series can be described in the 
real-frequency representation, and S(w) can be regarded 
as a functional of Gq(u>). Eqs. (27) and (28) relate the 
nonequilibrium quantities with the equilibrium ones. Es- 
pecially, at T = the usual zero-temperature formalism 
of the Green's function is applicable, and the right-hand 
side of eqs. (27) and (28) can be rewritten in terms of 
the vertex corrections. 14 '' Along this line, the low-energy 
behavior of the self-energy S r (w) and spectral function 
A(uo) = — ImG r (w)/7r has been calculated exactly up to 
terms of order oj 2 , T 2 , and (eV) 2 . 

4-2 Experimental determination of the Wilson ratio 

The low-energy behavior of the differential conduc- 
tance dl/dV has also been calculated up to terms of 
order T 2 and (eV) 2 using the results of A(ui) and the 
formula for the nonequilibrium current eq. (C-5) given 
in Appendix C. Specifically, in the electron-hole sym- 
metric case, the local Fermi liquid is characterized by 
two parameters, i.e., the Wilson ratio R and the energy 
scale A which corresponds to the width of the Kondo 
resonance. 23 ' In terms of these two parameters, the low- 
energy behavior of A{ui) and dl/dV can be expressed in 
the form 



7rA 



1-1 



{R-1) 



dl 
dV 



(R-lf 



h 



~a) 



3(R- 1) 



8 



eV 

X 



(30) 



l + 2(i?-l) 2 



TFT 

X 



5(i?-l) 2 feV 
A 



(31) 



These parameters contain all contributions of the pertur- 
bation series in U, and are defined with respect to the 
equilibrium ground state: R = Xs/l and A = A/7 where 
7 and Xs are the enhancement factors for the T-linear 
specific heat and spin susceptibility, respectively. 24 ) The 
values of these parameters can be evaluated using the 
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exact Bethe ansatz solution, 20, 25 > 26 ) as shown in Fig. 4. 
The width of the Kondo peak A decreases monotonically 
with increasing U and tends to A — > 4Tk/k for large U, 
where 



T K = ttA x /u/(2tt) exp[-7r 2 u/8 + l/(2u)] 



(32) 



and u — U I (ttA). The Wilson ratio increases with u from 
the nonintcracting value R = 1 , and converges rapidly for 
u > 2.0 to the value of the Kondo limit R — > 2 reflecting 
the suppression of the charge fluctuations by the strong 
Coulomb interaction. For the quantum dots, the value 
of R and A can be estimated experimentally from the 
results of dl/dV without carrying out the measurements 
of the spin susceptibility and specific heat: the coefficient 
of the T 2 and V 2 terms of dl/dV in eq. (31) may be 
estimated from the observations near the unitarity limit. 
This is also one typical feature of the Kondo system in 
the quantum dots. 

5. Summary 

Using the functional method, the dependence of the 
Keldysh Green's function G{lo) on eV and T has been 
confirmed to arise through the nonequilibrium distribu- 
tion function / e ff(w) which enters Gq(lu). From this 
property, the asymptotic behavior of G(u>) at high and 
low bias voltages has been deduced exactly. The low- 
voltage behavior is determined by a set parameters of the 
local Fermi liquid such as the Wilson ratio R and Kondo 
energy scale A. The values of these parameters can be 
estimated experimentally from the T 2 and V 2 contribu- 
tions of the differential conductance near the unitarity 
limit. In the high- voltage limit eV — > oo the Green's 
function becomes identical to that of the T — > oo limit 
at equilibrium, when the couplings between the dot and 
two leads are symmetric Tl =Tr. These results suggest 
that the Coulomb interaction U can be treated adiabat- 
ically at both low (eV < T K ) and high (U < eV) bias 
voltages. 
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Appendix A: Time Evolution of Density Matrix 



The time evolution of the density matrix is described 



by 



d_ 

dt 



-i[H,p(t)} . 



(AT) 



The formal solution of this equation can be obtained in 
the interaction representation, 



P(t) 



p(t)e 



-iHit 



= U{t,to)p{to)U(t ,t) , 



(A-2) 



where U(t,to) is the time-evolution operator 



U(t,t ) = Texp 



t 

i / dt'H 2 (t') 

'to 



(A-3) 



and H 2 (t) = e iHlt H 2 e~ iHlt The relations among the 
Schrodinger Os, interaction 0(t), and Hcisenberg On(t) 
operators are 



0(t) = e [Hlt O s e- [Hlt , 

o H {t)=u(o,t)d(t)u{t,o) . 



(AT) 
(A-5) 



To describe the nonequilibrium state, the initial condi- 
tion is given at to - > — 00 in eq. (A-2) assuming that each 
of the isolated leads is in the thermal equilibrium 



e -0[H^- l j, L N L -ij, R N R ] 
e -P[H 1 -n L N L -ij. R N Ii ] 



(A-6) 



for A — L, R. Then the cx- 



p(-oo) = 

where N x = Y, k o 
pectation value of 0#(t) with respect to p(0) can be 
written in the form of a time-ordered function along the 
Keldysh contour, 

(0 H (t)) = Tv[p(0)O H (t)} 

= Tr p(-oo)U(-cc,+oo)U(+oo,t)d(t)U{t,-oo) 

(A-7) 

Note that U(—oo, +oo) = W (+oo, — oo) can be expressed 
in terms of the anti-time-ordering operator T as 



U{— oo, +oo) = T exp 



J — ( 



dt'H 2 (t') 



(A-8) 



Appendix B: Matrix Dyson Equation 

In the matrix notation used in eq. (13), the Dyson 
equation for the intra-site Green's function at the dot 
G(to) and that for the inter-site ones between the dot 
and leads are given by 



G(lu) = g(uj)+g(Lj) S(w) G{w) 
+ v L g(uj)T 3 G Ld (uj) 

G xd (u;) = v x g x {bj)T 5 G{uj) , 
G dx (uj) = v x G(cj)T 3 g x (cj) , 



v R g(aj)r 3 G Rd (u>) , 
(B-l) 

(B-2) 

(B-3) 



T 3 = 



1 

-1 



P = T2 



1 

-1 



(B-4) 



Here the elements of G xd (t) and G xd (t) are defined 
by G Xd -(t) = -i(Tc Aff (t)4(0)), and G^(t) = 
— i (Td a (t) c\ a (0)), . . ., respectively. The Green's func- 
tions for the isolated leads g x (uj) include the distribution 
functions for the initial state, 



ffAM = p 

9 r x(<o) = E 





\<f> k x(r X )\ 2 
w - e kX + i0+ 



(B-5) 
(B-6) 
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ee [1 - 2 /^)MM -£>)]. (B-7) 

For the isolated dot, the Green's function can be writ- 
ten as {g(uj)}^ 1 = (uj — Ed) T3. Note that the contri- 
bution of the element g K (uj) which corresponds to eq. 
(B-7) vanishes for the isolated dot because of a property 
{g r (u)}- 1 g K (u){g a (u)}- 1 = 0. 11 ' Using eqs. (B-l) and 
(B-2), we have 

{G^)}" 1 = {GoMr 1 - H{u>) , (B-8) 

{G0M}- 1 = {ff(c)}- 1 - <r(u;) , (B-9) 

cr(w) = v 2 L T 3 g L (uj)T3 + v 2 R T 3 g R (uj)T 3 . (B-10) 

Then, using eqs. (B-5)-(B-7), <x(cj) is written in the form 



a a (uj) 



cr r (w) = vlg r L (uj)+v R g r R (uj) , 
n (w) - v\gf {w) + v 2 R g% (w) . 



(B-ll) 

(B-12) 
(B-13) 



When the density of states Img^(aj) is a constant and its 
band width is infinity, <r r (uj) becomes pure imaginary as 
v\ g\{uj) = — iT\. We obtain the explicit form of Go(w) 
given by eqs. (8)— (11) via eq. (B-9). Through eq. (B-8), 
G(uj) can be expressed in similar forms; 

G-(lo) = [1 - / off H]G>) + / off HG a M, (B-14) 
G-+(u) = - /cffH [G» - G» ] , (B-15) 
G+-( W ) = [1 - /cffH] [G r (uj) - G»] , (B-16) 
G ++ M = - [1 - / cff H]G» - f eS (cj) G», (B-17) 

where G r (ui) and / e s(w) are given by 
1 



oj- E d - a r (w) - S r (w) 



(B-18) 



~ /LHr £ + / fl Hr fl -^s-+H 
/offH " r £ + r fl -ims» ' (B49) 

Here f e s(u>) is real because S '"(w) is pure imaginary, 
and at equilibrium it becomes the usual Fermi function 
since E _+ (w)| v=0 = 2i/(w) Im E r (u>)| y=0 . n ) Note that 
S r = + There are some additional relations: 

£ a (w) = {S r (w)}*, and S~(w) = - {S++(w)}*. 



Appendix C: Current and Charge 

The operator for the current flowing from the left lead 
to the dot II and that from the dot to the right lead I R 
are given by 

I x = iew x J2 v x (4X-4c Aff ) , (C-l) 

<7 

where w R — 1 and wl = — 1 ■ The equation of continuity 
is written as dnd/dt + Ir — II =0 with rid = J2a n d<?- 



The averages of these operators with respect to p(0) are 



(n d ) = 2 y_ oo -G- +M 



(C-2) 



(7 A ) = 2ew A 



f 

j — ( 



do; 
2^ 



[G^M-G^M]- (C-3) 



Note that G^ and Gl^ can be expressed in terms of 
G using eqs. (B-2)-(B-3). Specifically, if the couplings 
with the leads satisfy a condition Tl{lu) cx T r (oj), the 
expectation value of the current can be rewritten as 27 ) 



r L (i R ) + r R (i L ) 
r R + r L 



(C-4) 



= ^ r d co[f L -f R ]^^[-lmG(c)]. (C-5) 
Note that (II) = (Ir) in the stationary state. 
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